A Modified Block Newton Iteration Forapproximating an Invariant Subspace of Asymmetric Matrixralf L

نویسنده

  • Gisela Timmermann
چکیده

In this paper we propose a Modiied Block Newton Method for approximating an invariant subspace S and the corresponding eigenvalues of a symmetric matrix A. The method generates a sequence of matrices Z (k) which span subspaces S k approximating S. The matrices Z (k) are calculated via a Newton step applied to a special formulation of the block eigenvalue problem for the matrix A, followed by a Rayleigh-Ritz step which also yields the corresponding eigenvalue approximations. We show that for suuciently good initial approximations the subspaces S k converge to S in the sense that sin ' k with ' k := ](S k ; S) converges to zero Q-quadratically provided that the eigenvalues belonging to S are separated from the remaining ones.

برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید

ثبت نام

اگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید

منابع مشابه

Certification Using Newton-Invariant Subspaces

For a square system of analytic equations, a Newton-invariant subspace is a set which contains the resulting point of a Newton iteration applied to each point in the subspace. For example, if the equations have real coefficients, then the set of real points form a Newtoninvariant subspace. Starting with any point for which Newton’s method quadratically converges to a solution, this article uses...

متن کامل

A New Inexact Inverse Subspace Iteration for Generalized Eigenvalue Problems

In this paper, we represent an inexact inverse subspace iteration method for computing a few eigenpairs of the generalized eigenvalue problem Ax = Bx [Q. Ye and P. Zhang, Inexact inverse subspace iteration for generalized eigenvalue problems, Linear Algebra and its Application, 434 (2011) 1697-1715 ]. In particular, the linear convergence property of the inverse subspace iteration is preserved.

متن کامل

A parallel Newton-Krylov flow solver for the Euler equations on multi-block grids

We present a parallel Newton-Krylov algorithm for solving the three-dimensional Euler equations on multi-block structured meshes. The Euler equations are discretized on each block independently using second-order accurate summation-by-parts operators and scalar numerical dissipation. Boundary conditions are imposed and block interfaces are coupled using simultaneous approximation terms (SATs). ...

متن کامل

Block Arnoldi-based methods for large scale discrete-time algebraic Riccati equations

Algebraic Riccati equations (discrete-time or continuous-time) play a fundamental role in many problems in control theory. They arise in linear-quadratic regulator problems, H ∞ or H 2-control, model reduction problems and many others. In this talk we propose numerical methods for large discrete-time algebraic Riccati equations (DARE). We present block projection methods that allow us to comput...

متن کامل

A Three-Dimensional Multi-Block Newton-Krylov Flow Solver for the Euler Equations

A three-dimensional multi-block Newton-Krylov flow solver for the Euler equations has been developed for steady aerodynamic flows. The solution is computed through a Jacobian-free inexact-Newton method with an approximate-Newton method for startup. The linear system at each outer iteration is solved using a Generalized Minimal Residual (GMRES) Krylov subspace algorithm. An incomplete lower/uppe...

متن کامل

ذخیره در منابع من


  با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید

برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید

ثبت نام

اگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید

عنوان ژورنال:

دوره   شماره 

صفحات  -

تاریخ انتشار 1996